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We study the deconfinement phase transition of compact U (1) pure lattice gauge theory with the Wilson action 
on closed topology lattices. In contrast to studies of compact QED on hypercubic lattices with periodic boundary 
conditions, we find no metastability signal at the phase transition on the lattices with the topology of a sphere. 
Thus the determination of the order of this phase transition has to be reconsidered. We argue that different 
properties of closed monopole loops on these topological inequivalent lattices might be responsible for the effect. 



1. INTRODUCTION 

For some time now computer studies of the 
phase transition of compact U(l) pure gauge 
theory with the Wilson action have produced first 
order signals. Original work [1] on comparatively 
small lattices could not resolve the energy gap, 
but subsequent high statistics Monte Carlo re- 
sults exhibited a metastability signal [2,3]. De- 
tailed studies of the latent heat in an extended 
parameter space corresponding to the adjoint re- 
presentation of the plaquette pointed towards the 
existence of a tricritical point at negative values of 
the adjoint coupling [4]. It was then generally ac- 
cepted that the compact U(l) pure gauge theory 
with the Wilson action has a weak first order PT 
with a small latent heat A(cos((7p)) ~ 0.016 in 
the thermodynamic limit. 

It was observed quite early [5] that this PT is 
accompanied by an intriguing behaviour of topo- 
logical excitations, so-called monopoles. These 
are artifacts of the compact formulation of the 
(7(l)-action, i.e. fiux deficits in 3-cubes. On the 
dual lattice they form closed loops due to cur- 
rent conservation. In fact, suppressing them af- 
fects the PT, moving it to the disordered confi- 
nement phase [6] . Moreover eliminating monopo- 
les from the path integral by introducing certain 
constraints makes the confinement phase disap- 
pear for positive values of /3 [7] . 

Almost all numerical simulations have used pe- 
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riodic boundary conditions on hypercubic latti- 
ces, some have chosen helical b.c, which do not 
behave differently with regard to the subject dis- 
cussed here. Hypercubic lattices with periodic 
boundary conditions are homotopic the surface 
of a torus with the fundamental homotopy group 
TTi ~Z®Z®Z®Z. The formulation of compact 
QED on the torus contains certain peculiarities 
with respect to the monopole excitations of the 
theory: A monopole loop (on the dual lattice) 
may be considered as the boundary of a Dirac 
sheet. Whereas the shape of the Dirac sheet is 
gauge variant its boundary is not. A pair of mo- 
nopole loops wrapping around the torus may be 
connected by a Dirac sheet and pairwise anni- 
hilation may leave a closed Dirac sheet in the lat- 
tice configuration. This has been observed indeed 
on finite lattices as a metastable energy gap in the 
Coulomb phase of the model [8] . Furthermore due 
to the periodic b.c. one finds monopole loops that 
are closed only because of the periodic boundary 
conditions and that may not be continuously de- 
formed (are not homotopic to) to points. They 
are present in the "hot" phase of the model, but 
usually not in the "cold" phase. These particular 
properties of the theory on a torus might have 
led to a wrong identification of the order of the 
QED deconfinement PT with the Wilson action. 
In fact earlier MCRG studies of the model sho- 
wed in the vicinity of the Wilson line a fiow of 
couplings, which resembled the fiow diagram of a 
second order critical point quite well [9] . 
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2. LATTICE WITH TRIVIAL HOMO- 
TOPY GROUP 

The central idea is to study compact QED with 
the Wilson action on a 4D lattice, which has the 
topology of a sphere 7ri(5'4) ~ 1. On such a lat- 
tice any closed monopole loop is homotopic to a 
point and therefore above mentioned complicati- 
ons due to the periodic boundary conditions are 
excluded. For the construction we choose the D- 
manifold as the boundary of a _D-|- 1-manifold, for 
which we select a D + 1-dimensional hypercube 
Ho+i{L) with linear extend L and L^'^^ sites. 
We denote its boundary by SHo (L); it consists of 
all sites which have at least one cartesian coordi- 
nate value 1 or _L. In Table 1 we display the num- 
ber of basic elements (sites, links, plaquettes) of 
the manifold 5'-ff4(i.). Ratios of e.g. ^links/'^sites 
approach for the SH4(L) lattice for large L asym- 
ptotically the values for the hypercubic geometry 
H4L). Whereas Hd(L) is self-dual, SHd(L) is 
obviously not. It is however possible to construct 
the dual manifold SH'j^ (L) (which has the same 
homotopy properties) and we have used this con- 
struction for some exploratory identifications of 
closed monopole loops. 



Table 1 

Characteristic numbers on SH4(L) lattices. 



manifold 




SH4{L) 


sites 
links 
plaquettes 


10(1. - 
40(1. 
60(1. 


1)4 + 20(1.-1)2+2 

- 1)4 + 40(1.-1)2 

- 1)4 + 20(1.-1)2 



3. SIMULATION AND RESULTS 

We have simulated compact QED with the Wil- 
son action on lattices SH (A) , SH (Q) , SH (S) and 
SH(IO). Comparing the number of degrees of 
freedom, i.e. the number of link variables, we 
find that those systems roughly correspond to hy- 
percubic lattices _ff(6), -ff(9), -ff(12) and _ff(16) in 
D = 4. The geometrical properties have been im- 
plemented by tables and we use a 3-hit Metropolis 
update. 



We have studied long runs of (9(10®) sweeps 
at various couplings /3 close to the deconfine- 
ment PT. We then combine our data using the 
Ferrenberg-Swendsen technique and obtain an op- 
timal estimator for the density of states. From 
that we determine various quantities like the in- 
ternal energy, the specific heat Cv{P,L) and the 
BCL-energy-cumulant V{p,L). 

Naively one expects any differences due to the 
topology of the lattice to disappear in the large 
volume limit. Therefore we expect similar pseu- 
docritical coupling estimators for our case than 
in the case of the torus topology. The pseudocri- 
tical couplings as defined by the position of the 
maxima of the specific heat as a function of L are 
/3e(4) ^ 1.0155, /3c(6) ^ 1.0135, /3c(8) ^ 1.0127 
and /3c(10) ~ 1.0122. Simulations for the theory 
on the hypercubic lattice with periodic boundary 
conditions show clear two-state signals with a pro- 
nounced double peak distribution in the energy 
near the PT. States with energies in between 
the pure phases are clearly suppressed. In our 
runs, however, we could not find any such signals. 
Away from the pseudocritical points the energy 
distribution histograms always were purely Gaus- 
sian, while at pseudocriticality their shape some- 
what departs form a Gaussian form. Fig. 1 com- 
pares the energy distribution functions of a run 
on the SH(S) lattice with a run on a -ff(lO) lat- 
tice. In both runs the couplings were chosen very 
close to the peak locations of the specific heat. 

We did however observe very long autocorre- 
lation times up to (9(1000) sweeps for our lar- 
gest lattices at pseudocriticality. We also analy- 
zed the finite size scaling properties of the pseu- 
docritical couplings (3c{L) as well as the scaling 
properties of the maximum of the specific heat 
Cy^'^iL) and the minimum of the BCL-energy- 
cumulant V™"{L). Fitting these quantities in 
the i.-interval 6 < i. < 10 with the large L sca- 
ling laws 

we obtain the following rough estimates for — 
3.1, V K l.<d , y™'" ^ 0.66655 < | and ^ 
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Figure 1. A comparison of energy distribution hi- 
stograms for compact QED with the Wilson ac- 
tion for the -ff(lO) hypercubic lattice with peri- 
odic boundary conditions and the SH(S) closed 
topology lattice. The simulation on -ff(lO) have 
been done in the multicanonical Ensemble. The 
accumulated statistics on the SH(S) lattice is lar- 
ger than half a million sweeps. 



1.0102. These estimates should be taken with 
precaution. It might well be that the asymptotic 
scaling behaviour for above finite size scaling laws 
sets in only on very large lattice sizes. While in 
general not much is known about this for gauge 
theories, simpler models, like e.g. D = 2 Potts 
models require lattices as large as 10 times the 
finite correlation length at the first order PT in 
order to show the above asymptotic scaling beha- 
viour [10]. It is therefore no surprise that those 
estimates are contradictory among themselves: a 
value of I indicates a first order PT, 

while the value of the Cy-exponent ^ ^ 4 does 
not. 

Right now our conclusion has to be that we 
have not yet sufficient statistics and large enough 
lattice sizes to unambiguously decide on the order 
of the PT from the simulations for the Wilson ac- 



tion. However the blunt absence of metastability 
signal in our data, as compared to the standard 
case, makes it hard to accept that a first order 
signal will appear in the thermodynamic limit. 
After all, first order phase transitions are defined 
by coexistence of metastable bulk states forming 
an interface in between them. Our data show no 
sign of such an interface (cf. Fig. 1). As obvi- 
ously now there is reason to doubt the usual iden- 
tification as first order PT we think, that due to 
the interest on other phase transitions of theories 
with the compact (7(l)-gauge group (like QED 
with matter fields) this problem should be clari- 
fied. Like emphasized in the studies which pre- 
dicted the tricritical point of QED [4], a possible 
approach to the resolution of the problem would 
be the inclusion of adjoint plaquette terms into 
the action. Simulations in this direction are on 
the way. 
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